import scipy.stats as stats
import numpy as np

mu = 4 #5分钟有电话次数
k = 3 #接电话3次
X = np.arange(0, k+1, 1)
plist = stats.poisson.pmf(3, mu)
print('5分钟内接到3次电话的概率:',plist[3]) #5分钟内接到3次电话